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Inflationary cosmology leads to the picture of a "multiverse," involving an infinite number of 
(spatially infinite) post-inflationary thermalized regions, called pocket universes. In the context of 
theories with many vacua, such as the landscape of string theory, the effective constants of Nature are 
randomized by quantum processes during inflation. We discuss an analytic estimate for the volume 
distribution of the constants within each pocket universe. This is based on the conjecture that the 
field distribution is approximately ergodic in the diffusion regime, when the dynamics of the fields is 
dominated by quantum fluctuations (rather than by the classical drift). We then propose a method 
for determining the relative abundances of different types of pocket universes. Both ingredients are 
combined into an expression for the distribution of the constants in pocket universes of all types. 



I. INTRODUCTION 

The fundamental theory of nature may admit multiple 
vacua with different low-energy constants. If there were 
just a few vacua, as in standard GUT models, then a few 
observations would determine which one corresponds to 
the real world. Predictions would then follow for every 
other observable in the low energy theory. However, it 
has recently been realized that in the context of string 
theory there may be a vast landscape of possibilities, with 
googols of vacua to scan 0,0, HI ■ Many of these may look 
very much like our own vacuum, except for slight vari- 
ations in the values of the constants. On the surface, 
this seems to undermine our ability for predicting these 
values, even after a systematic examination of the land- 
scape. 

Cosmology, on the other hand, suggests that rather 
than giving up on our ability to make predictions, we 
may in fact broaden their scope. Thus, instead of trying 
to determine from observations which vacuum is ours, 
we may try to determine, from the theory, what is the 
probability for the observation of certain values of the 
constants. 1 Indeed, eternal inflation 0, [j| leads to the 
picture of a "multiverse," where constants of nature take 
different values in different post-inflationary regions of 
spacetime. Observers bloom in such thermalized regions 
at places where the conditions are favorable, much like 
wildflowers at certain spots in the forest. Given a refer- 
ence class of observers, we can ask what is the probability 
distribution for the values of the constants that they will 
measure. This approach was suggested in and further 



1 Note that this question is relevant even if the number of vacua 
is small. 



developed in 0,11,11 03 It leads to the following formal 
expression for the probability of observations, 

Pobs(X) ozP(X) n ohs {X). (1) 

P(X) is the volume fraction 2 of thermalized regions with 
given values of the constants X , and n b s (X) is the num- 
ber of observers in such regions per unit thermalized vol- 
ume. 3 

Even though the dynamics responsible for the random- 
ization of the constants during inflation is well under- 
stood, the calculation of probabilities has proven to be 
a rather challenging problem. The root of the difficulty 
is that the volume of thermalized regions with any given 
values of the constants is infinite (even for a region of 
a finite comoving size). To compare such infinite vol- 
umes, one has to introduce some sort of a cutoff. For 
example, one could include only regions that thermal- 
ized prior to some time t c and evaluate volume ratios in 
the limit t c — > oo. However, one finds that the results are 
highly sensitive to the choice of the cutoff procedure (in 
the example above, to the choice of the time coordinate 



2 P(X) is sometimes referred to as the "prior distribution." This 
name is motivated by Carter's original discussion of anthropic se- 
lection in terms of Bayes' rule. To avoid confusion with the usage 
of prior distributions in other contexts, here we shall simply call 
it the thermalized volume distribution. Note that the concept 
of a prior distribution for a set of parameters is also used when 
fitting observational data to a given model. There, the "prior" is 
often no more than a guess, representing our ignorance or prej- 
udice about, say, the allowed range of the parameters. On the 
contrary, P(X) is here a quantity which should be calculable 
from the theory. 

3 n obs{X) is often referred to as the anthropic factor. Readers who 
dislike anthropic arguments are encouraged to think in terms of 
reference classes of their own choice. 
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FIG. 1: Structure of eternally inflating spacetime. Thermal- 
ized regions are shown by dark shading. 



t 0,0]; see also 0,^1 for a recent discussion.) The 
reason for the cutoff dependence of probabilities is that 
the volume of an eternally inflating universe is growing 
exponentially with time. The volumes of regions with 
all possible values of the constants are growing exponen- 
tially as well. At any time, a substantial part of the total 
thermalizcd volume is "new" and thus close to the cutoff 
surface. It is not surprising, therefore, that the result 
depends on how that surface is drawn. 

As suggested in 0, @, this difficulty can be circum- 
vented by switching from a global distribution, defined 
with the aid of some global time coordinate, to a dis- 
tribution based on individual thermalized regions. The 
spacetime structure of an eternally inflating universe is 
illustrated in Fig. ^ The vertical axis is the proper 
time measured by comoving observers, and the horizontal 
axis is the comoving coordinate. Thermalized regions are 
marked by grey shading. Horizontal slices through this 
spacetime give "snapshots" of a comoving volume at dif- 
ferent moments of (global) time. Initially, the whole vol- 
ume is in the inflating state. While the volume expands 
exponentially, new thermalized regions are constantly be- 
ing formed. These regions expand into the inflating back- 
ground, but the gaps between them also expand, making 
room for more thermalized regions to form. The thermal- 
ization surfaces at the boundaries between inflating and 
thermalized spacetime regions are 3-dimensional, infinite, 
spacelike hypersurfaces. The spacetime geometry of an 
individual thermalized region is most naturally described 
by choosing the corresponding thermalization surface as 
the origin of time. The thermalized region then appears 
as a self-contained infinite universe, with the thermaliza- 
tion surface playing the role of the big bang. Following 
Alan Guth, we shall call such infinite domains "pocket 
universes." All pocket universes are spacelike-separated 
and thus causally disconnected from one another. In 
models where false vacuum decays through bubble nucle- 
ation, the role of pocket universes is played by individual 
bubbles. 



The proposal of Refs. 0, II] applies to the special 
case when all pocket universes are statistically equiva- 
lent. This happens, for instance, if there is only one 
low-energy vacuum, and the variation of the constants 
X is due to long-wavelength fluctuations of some nearly 
massless fields. In this case, one may use any one of 
these pockets in order to calculate P(X), which is then 
proportional to the volume fraction occupied by the cor- 
responding regions in the pocket universe. This fraction 
can be found by first evaluating it within a sphere of large 
radius R and then taking the limit R — > oo. 

Although the above proposal is conceptually satisfac- 
tory, it remains unclear how it should be implemented in 
practice. The most direct method, a numerical simula- 
tion, runs into severe computational limitations 0, 1 1 1| , 
and analytic methods have been developed only for very 
special cases. Moreover, if two or more vacua are mutu- 
ally separated by inflating domain walls, then they can- 
not coexist in the same pocket universe. In this situation 
(which is expected to be quite generic in the landscape 
of string theory) there are distinct types of pocket uni- 
verses, and we have to face the problem of comparing 
probabilities of different pockets. The purpose of the 
present paper is to try to improve on the approach of 
Refs. 7, 8, 9] by addressing both of these issues. 

In Section II we discuss an analytic estimate for the 
volume distribution of the constants within each pocket, 
based on the conjecture that the field distribution is ap- 
proximately ergodic in the diffusion regime. When there 
are different types of pockets, the same approach can 
be used in order to find the internal volume distribution 
P(X; j) in a pocket of type j. Then, in Section III, we in- 
troduce a new object pj which characterizes the relative 
abundance of each type of pocket universe. The aficiona- 
dos may be aware of previous exploratory definitions of 
Pj, which were nevertheless afflicted with certain draw- 
backs 10]. For comparison, in Section IV we discuss one 
of these alternatives and in Section V we illustrate both 
definitions with some examples. Section VI discusses how 
the two objects P(X;j) and pj may be combined in or- 
der to find the full distribution for the constants. Our 
conclusions are summarized in Section VII. 



II. PROBABILITIES WITHIN A POCKET 
UNIVERSE 

A. Symmetric inflaton potential 

The spacetime structure of a pocket universe is illus- 
trated in Fig. [5] The surface X* in the figure is the ther- 
malization surface. It is the boundary between inflating 
and thermalized domains of spacetime, which marks the 
end of inflation and plays the role of the big bang in 
the pocket universe. The surface S g is the boundary be- 
tween the stochastic domain, where the dynamics of the 
inflaton field is dominated by quantum diffusion, and the 
deterministic domain, where the dynamics is dominated 
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quantum diffusion 



FIG. 2: Structure of a pocket universe. 



by the deterministic slow roll. Thus, E g marks the onset 
of the slow roll. 

The surface £* is always spacelike. We shall assume 
that T, q is spacelike as well. This can be arranged if 
T, q is chosen not right at the diffusion boundary, but 
somewhat into the slow roll domain. We shall comment 
on the precise definition of £ g at the end of Section II. C. 

The probability P(X)dX for a variable X is defined as 
the volume fraction on the thermalization hypersurface 
£* with given values of X. It can be expressed as 



P(X) <x P q (X)Z 3 (X), 



(2) 



where P q (X) is the distribution (volume fraction) of X 
on E g , and Z(X) is the volume expansion factor during 
the slow roll in regions with the value X. For simplicity, 
we shall identify the variable X with the scalar field re- 
sponsible for its value. Here, we have assumed that the 
diffusion of X can be neglected during the slow roll. We 
shall assume also that X interacts very weakly with the 
inflaton, so that it does not change appreciably during 
the slow roll. Otherwise, we would have to include an 
additional Jacobian factor 



det [dX(t q )/dX(t obs 



(3) 



in Eq. J3J) g|- Here, t q and t b s indicate, respectively, the 
time when slow roll begins and the time when observa- 
tions are made. 

The expansion factor Z(X) can be found from 



Z{X) w exp 



4-7T 



y " x H(<p,X) 



d(t> 



where cj> is the inflaton field, 

H(c/>;X) = [8wV(cl>,X)/3] 1 / 2 



(4) 



(5) 



is the inflationary expansion rate, V(4>, X) is the inflaton 
potential, the prime stands for a derivative with respect 
to (f>, and we use Planck units throughout this paper. We 
denoted by <fi q x and 4>*x the values of 4> at the bound- 
ary surfaces Ti q and £*, respectively. These values are 
defined by the conditions 



and 



H'/H 2 {0 qX ,X) ~ 1 



H'/H{<i>*x,X)~l. 



(6) 



(7) 



The subscript X in <p q x indicates that the value of <p 
where slow roll begins is generally A-dependent, <f> = 
<f>q(X) = <f> q x, and similarly for <p*x- 

The distribution P q (X) can in principle be determined 
from numerical simulations of the quantum diffusion 
regime. Some useful techniques for this type of simu- 
lation have been developed in However, the 
simulations quickly run into computational limits, due 
to the exponential character of the expansion. Analytic 
techniques for the calculation of P(X) are therefore very 
desirable. 

A special case where the analysis is trivial is the class 
of models where the potential is independent of X in the 
diffusion regime. Then P q (X) can be determined from 
symmetry, 



P q (X) = const. 



(8) 



An example is a "new" inflation type model with a com- 
plex inflaton, <p = \4>\ exp(iA). Inflation occurs near the 
maximum of the potential at <f> = 0, and we assume in 
addition that the potential is symmetric near the top, 
V = V(\(j)\). Equation (JSJ follows if this property holds 
throughout the diffusion regime. In such models, the dis- 
tribution (0 reduces to 



P(X) cx Z 3 {X). 



(9) 



Thus, the volume distribution is simply determined by 
the expansion factor during the slow roll. 

In models where the inflaton potential does not have 
the assumed symmetry, the factor P q (X) may provide 
an additional dependence on X. We now turn to the 
discussion of this more general case. 



B. Eternal inflation without thermalization 

Let us first consider a model of eternal inflation with- 
out any low-energy vacua, where the dynamics is dom- 
inated by diffusion everywhere in the potential. The 
stochastic evolution of fields in the course of eternal in- 
flation is described by the distribution function 



F(4>a,t), 



(10) 
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where <p a stands collectively for the inflaton and other 
light fields. This distribution satisfies the Fokker- Planck 
(FP) equation 



where 



dF = dJg 

dt d(j> a ' 
where the flux J a (4>, t) is given by 

J a = -^-[D{^)F]+v a {cp)F. 



Here. 



D(<(>) = H((t>f +2 /8Tr 2 



is the diffusion coefficient, 



(11) 



(12) 



(13) 



(14) 



is the "drift" velocity of the slow roll, and as before 
H(<f>) = [87rT^((/))/3] 1 / 2 is the inflationary expansion rate. 

The parameter (3 in Eqs. ((13(1 and (|f 4|l represents the 
freedom of choosing the time variable t, which is assumed 
to be related to the proper time of comoving observers r 
by 



dt = [H MrW-Pdr, 



(15) 



where the functions cj>(r) are taken along the observer's 
world lines. Hence, (3 — I corresponds to the proper time 
parametrization t = r and (3 = corresponds to using the 
logarithm of the scale factor as the time variable. (We 
shall refer to the latter choice as the scale factor time.) 

The FP equation is usually supplemented by boundary 
conditions at the thermalization boundary in </>-space, 
S*, where the conditions of slow roll are violated, and in 
some cases at the Planck boundary S p , where V(tf>) ~ 1, 
and quantum gravity effects become important. Here, we 
shall assume that neither of these boundaries is present. 
For example, the </>-field space could be compact, with 
the potential V(4>) satisfying V{<j>) < f and H'/H 2 < f 
everywhere in this space. Then the fields <p a will drift 
constantly from one value to another. The corresponding 
distribution has a stationary form with a vanishing flux, 



Ja =0. 



(16) 



With J a from iff^l . the solution of Eq. (fTB|> is easily 
found Up to a normalization constant, it is given 

by 



(17) 



Let us now use the scale factor time, (3 = 0, which cor- 
responds to measuring time in units of the Hubble time 
if -1 . With this choice, Eq. 1(17(1 can be rewritten as 



S(</>) = it/ H 2 (cj)) 



(19) 



is the Gibbons-Hawking entropy of de Sitter space. 

The solution (|f 8(1 has a simple physical interpretation. 
The distribution F((j>, t) can be thought of as represent- 
ing observations by a comoving observer. The observer 
can see only her own horizon region (/i-region), and the 
value of <p should be understood as an average over that 
region. The quantum state of the /i-region is constantly 
changing due to de Sitter quantum fluctuations, and the 
average <j> is also changing as a result. The distribution 
function F((f>) is proportional to the fraction of time spent 
in quantum states with a given <fi. The number of micro- 
scopic quantum states corresponding to a given value of 
(j> is exp[S(<f)], and Eq. (|f 8fl indicates that, apart from 
the prefactor, the evolution during eternal inflation is er- 
godic, with all quantum states explored with an equal 
weight. 

The prefactor H~ 2 suggests that at higher rates of in- 
flationary expansion, the quantum state of an /i-region 
changes at a higher pace. The scale-factor time spent 
in each quantum state is St oc H~ 2 (the proper time 
St oc H~ 3 ). We note, however, that the prefactor in the 
solution (|f 8(1 depends on the choice of ordering of the 
non-commuting factors <f> and d/d(j) in the FP equation. 
In Eq. 1(12(1 above, we assumed the Ito factor ordering — 
the choice suggested by the phenomenological analysis in 
[13. A full microscopic derivation would require inclu- 
sion of quantum gravitational fluctuations, and until then 
the precise form of the prefactor will remain uncertain. 
We shall mostly disregard the prefactor in what follows. 

The distribution 1(18(1 has an alternative interpretation, 
which will also be important for our analysis. It follows 
from the so-called physical- volume form of the FP equa- 
tion 111 111 111, 



dt d(pa 



(20) 



Here, the function F((f>, t) characterizes the distribution 
of physical volume between regions with different values 
of <j>. More precisely, F(<j>, t)d n (f> is the fraction of volume 
occupied by regions with <fi a in the intervals d<p a on hy- 
persurfaces of constant t. The flux J a is given by 1(12(1 
with F replaced by F, and the last term in l(2U|) accounts 
for the expansion of physical volume. 

With = 0, solutions of Eq. and those of Eq. 
are closely related > 



F(cf>,t)=e M F(<p,t) 



(21) 



F^) = H- 2 (^)e s ^\ 



(18) 



This applies, in particular, to the solution 1(18(1 . On sur- 
faces of constant t, the factor exp(3£) is a constant, and 
thus the volume distribution of cf> along these surfaces is 
the same as the time distribution along the worldlines of 
comoving observers. 

To appreciate the significance of this result, imagine 
dividing the spacetime into AD cells of size ~ S/H with 
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6 -C 1. This can be done by first foliating the spacetime 
by surfaces of constant scale factor time, t — 0,(5,2(5, .... 
The timelike separation between consecutive surfaces is 
of order 6/H. The next step is to divide the surface 
t = into 3D cubes and to extend these cubes along 
comoving geodesies all the way to the surface t = 6. As 
a result, the spacetime layer between the surfaces t = 
and t = 8 is divided into AD cells. This procedure is 
similar to laying bricks, starting with the surface t = 0, 
except that the bricks are somewhat irregular, with their 
tops being slightly bigger than their bottoms and their 
size depending on the local value of H. 

Now, we wish to know how cj> is distributed among dif- 
ferent cells. It follows from Eq. I|21|) that the distribution 
is given by l|18fl. both in the time direction and in the 
spacelike directions along the surfaces t = const. We can 
expect, therefore, that the 4-dimensional distribution of 
4> among the cells, that is, the distribution in a randomly 
picked AD volume, is still given by (|18fl . Note that the 
form of the distribution does not depend on the choice 
of the spacelike hypersurface t — 0. This indicates that 
the distribution is not sensitive to the details of how the 
spacetime is divided into cells, as long as the cell sizes 
are set by the local horizon. 



C. Ergodic conjecture 

We now turn to the more physically interesting case, 
when inflation can end by thermalization in a low-energy 
vacuum. Here we shall assume that there is only one 
such vacuum, so that all thermalized regions are statis- 
tically equivalent. Each thermalized region has its own 
infinite hypersurface £ 9 , which marks the boundary be- 
tween quantum diffusion and slow roll. We are interested 
in the distribution P q (X) on E g . Here, we have switched 
back from the notation <j> a to (</>, X), where X stands for 
one or several light fields. 

The simplest guess is that P q (X) has the same form as 
the distribution (|18fl in the diffusion region, 

P q {X)^H- 2 {4> qX ,X)eMS{4> q x,X)]. (22) 

Here, as before, (f> q x is the value of <f> a t the onset of the 
slow roll for given values of the fields X. We do not have 
a proof that this guess is correct in general, but there are 
some special cases where the distribution (1221) does seem 
to apply. 

Suppose that the inflaton 4> is trapped in a metastable 
false vacuum, <f> — <pp- The false vacuum decays through 
bubble nucleation, which is followed by slow roll and ther- 
malization, as in models of "open inflation" [23]. Prior 
to bubble nucleation, the evolution of the fields X is gov- 
erned by the stochastic dynamics described in the previ- 
ous subsection. The corresponding probability distribu- 
tion is 

P(X) cx H- 2 ((/) F ,X)exp[S((t) F ,X)]. (23) 



Once a bubble nucleates, the bubble wall expands 
rapidly with a speed approaching the speed of light, so 
the worldsheet of the wall closely follows the future light- 
cone of the nucleation point. The onset of slow roll sur- 
face Tj q will follow the same lightcone on the interior side 
of the bubble. If we imagine the spacetime divided into 
horizon-size cells with different values of X, we can ex- 
pect that the bubble will cut through a representative 
sample of cells and that the distribution of X on S 9 will 
be given by (|23|l . The distribution on the thermalization 
surface £* is obtained by multiplying i|23|) with the slow 
roll volume expansion factor, as in Eq. @. 

It should be understood that if the fields X interact 
with the bubble, then the distribution (|2*3*|l is processed as 
it goes through the bubble walls. This effect should be in- 
corporated, just as the evolution of the fields during slow 
roll is incorporated in the Jacobian . In certain cases, 
the result of this processing can be easily estimated. For 
fields X whose mass is much smaller than the inverse size 
of the bubble at the time of nucleation, the distribution 
inside the bubble takes the form [2l| 

P q {X) ~ e- lB{x \ (24) 

Here, the subindex q is used to denote that this is the dis- 
tribution at the beginning of slow roll inside the bubble, 
right after nucleation. Ib{X) in the exponent is the ac- 
tion of the Coleman-de Luccia instanton describing bub- 
ble nucleation, and an adiabatic approximation is used 
for the Euclidean solution, where the field X is assumed 
to take a constant value. In fact, all values of X are real- 
ized at distant places inside any given bubble, but their 
distribution is still given by Eq. l(2"4T 21] . Note that in 
the limit where there is no bubble, the Euclidean action 
is minus the entropy of de Sitter space, Ig = —S, and 
the expression <|24[) reduces to (|23H . 

In cases where both <f> and X fields are in the diffusion 
regime prior to the slow roll, the situation is less clear. 
However, the following heuristic argument suggests that 
the distribution i|22fl may still approximately apply. 

Consider an ensemble of comoving observers who start 
their evolution in regions where the fields (0, X) are near 
the top of their potential. Most of these observers will 
get to the slow roll and thermalization after a relatively 
short period of diffusion. Only a small fraction of the ob- 
servers will stay in the diffusion regime for much longer, 
but their respective regions will be expanded by a huge 
factor. This is the essence of eternal inflation: the histo- 
ries of these "atypical" observers are, in a certain sense, 
more representative than those of the typical ones (af- 
ter all, the thermalization surface is infinite thanks to 
the contribution of regions which linger in the diffusion 
regime for an indefinite amount of time.) 

The measurements made by observers who spend a 
long time in the diffusion regime are described by the FP 
equation (|llf) . supplemented by the constraint that the 
boundary S q between the diffusion and slow roll regions 
in the field space cannot be crossed. [The surface S q is 
defined by the condition Mathematically, this can 



6 



be enforced by imposing the reflecting boundary condi- 
tion on S q . The solution of the FP equation is then given 
by the stationary distribution (JTBJl. This indicates that, 
during diffusion, the observers will see ergodic evolution 
in their respective /i-regions. 

We can thus argue, as in the previous subsection, that 
Eq. I|18|) describes the AD distribution in the quantum 
diffusion region of spacetime, which is in the past of the 
hypersurface E q . Then it appears plausible that the same 
distribution may extend, at least approximately, to the 
boundary S g with the slow roll region. This is the ergodic 
conjecture (|22|) . It assumes that all quantum states corre- 
sponding to the field values on S q are equally represented 
by /i-regions on S g . 

It should be noted that <f> q x is defined by Eq. 10 only 
in the order-of-magnitude sense. This results in a sig- 
nificant uncertainty 4 in the distribution l|22[l . due to the 
exponential dependence on S((f> q x, X). A possible way of 
defining S g more precisely is the following. Suppose first 
there is a single inflaton cf> and no other light fields. Sur- 
faces of constant (f> are spacelike in the slow roll regime, 
but not necessarily so in the diffusion regime. For <fi in the 
diffusion regime, the surfaces have a complicated topol- 
ogy, while in the slow roll regime the topology is trivial. 
The transition between the two might be sharp, having 
the character of a phase transition. We could then define 
4> q to be the critical value of <fi. If there are several fields, 
the phase transition on the surfaces H' /H 2 — const could 
be used to define S q . 

D. An example 

As an illustration, let us consider the potential 

V{^X) = \m\X)^\ (25) 

where m(X) is a very slowly varying function of the fields 
X. Then 

H{l>,X) = {-K/c)m{X)<t>, (26) 

where c = (37r/8) 1/2 , and Eq. Q gives 0* — I. 

To make the condition © more precise, we rewrite it 

as 

H'/H 2 (<j) qX ,X) = k (27) 
with k ~ 1, which yields 




We shall assume that m(X) -C I; then 4>* *C <fi q . 



4 This uncertainty is not present in models with bubble nucleation 
discussed earlier in this subsection. 



The ergodic factor in Eq. (|22|l is given by 

exp[7r/tf 2 (0 8 x, X)} = exp[ck/m(X)}, (29) 

and the volume expansion factor is 

Z 3 « exp(67r<^ x ) = exp[6e/fcm(X)]. (30) 

We see that both factors are sensitive to the value of k 
chosen to define <fi q x- On the other hand, both factors 
have the form exp[C/m(X)] with C ~ 1, so we can be 
reasonably confident that the full distribution has this 
form, with a strong peak at the smallest allowed value of 
m(X), 

III. COUNTING POCKETS 

In models with pockets of several different types, we 
need to consider an additional factor characterizing the 
relative abundances of different types of pockets. Con- 
sider, for instance, a false vacuum which can decay into 
true vacua 1 or 2, with respective tunneling rates K\ and 
«2- If «i > K2, it seems intuitively clear that there will 
be more pockets of type I than of type 2 in the multiverse. 
This notion, however, should be made more precise. 

In this section we propose a method for determining 
the relative abundances pj of different types of pockets. 
The counting is performed at the future boundary of 
spacetime — the place where everything has been said and 
done. Of course, the number of pockets of any given type 
is infinite, and a cutoff is needed. We shall argue that 
one such cutoff suggests itself naturally, and this allows 
us to perform an explicit calculation of the pj 's. One can 
then examine different examples to check whether this is 
a reasonable prescription. As mentioned in the introduc- 
tion, there have been in the past alternative proposals 
for pj which nevertheless are not very satisfactory. For 
comparison, one of these will be discussed in Section IV. 

The full probability distribution should be obtained by 
combining the weights pj with the volume distributions 
P(X;j) in a suitable way. This will be discussed in Sec- 
tion VI. 



A. Models with bubble nucleation 

The string theory landscape is expected to have a mul- 
titude of high-energy metastable "false" vacua. Such 
vacua can decay through bubble nucleation. Bubbles of 
lower-energy vacuum nucleate and expand in the high- 
energy vacuum background. If the "daughter" vacuum 
has a positive energy density, then inverse transitions are 
also possible: bubbles of high-energy vacuum can nucle- 
ate in the low-energy one [22] • (This is the so-called "re- 
cycling" process.) In both cases, the radius of the bub- 
bles asymptotically approaches the comoving horizon size 
in the parent vacuum at the moment of nucleation |23| . 
In general, we will have multiple bubbles within bubbles 
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within bubbles of many different types. The endpoints 
of this evolution are the negative or zero-energy vacua, 
which do not recycle. We shall call them terminal vacua. 

Let fj (t) be the fraction of comoving volume occupied 
by vacuum j at time t. The evolution equations for fj(t) 
can be written as 



dt 



(31) 



where the first term on the right-hand side accounts for 
loss of comoving volume due to bubbles of type i nucle- 
ating within those of type j, and the second term re- 
flects the increase of comoving volume due to nucleation 
of type-j bubbles within type-z bubbles. The transition 
rate Kij is defined as the probability per unit time for an 
observer who is currently in vacuum j to find herself in 
vacuum i. It can be expressed as 



(32) 



where Ajj is the bubble nucleation rate per unit physical 
spacetime volume and Hj is the expansion rate in vac- 
uum j. At this point, it will be convenient to distinguish 
between the inflating vacua, which we will label by Greek 
letters, and the terminal ones, for which we will reserve 
the indices m and n. Then, by definition, 



A, 



0. 



(33) 



We are interested in the number Nj of nucleated bub- 
bles of a given low energy vacuum j (not necessarily a 
terminal one). As we know, the number of bubbles of 
all kinds grows without bound, even within a region of a 
finite comoving size. We thus need to cut off our count. 
The method we propose is to include only bubbles greater 
than some small comoving size e. We will take the limit 
e — > in our final result. Note that this cutoff prescrip- 
tion is independent of time parametrization. The comov- 
ing bubble size is set by the horizon at the moment of 
nucleation and remains constant in time. Our bubble 
counting can therefore be performed at future infinity. 
We shall call this the comoving horizon cutoff (CHC) 
method. 

We can relate the number of bubbles nucleated in an 
infinitesimal time dt in the parent vacuum a, dNj a (t), to 
the accompanying increase in comoving volume dfj a = 
Kjafadt (which is the product of the number of bubbles 
and the comoving volume each bubble covers) : 



Ait 

dN ja (t)—H, 



- 3 a- 3 



(t) — Kj a f a dt. 



(34) 



Here, H a is the Hubble constant of the parent vacuum, 
a(t) is the scale factor, and we have used the fact that 
the comoving radius of bubbles of any type j nucleating 
in vacuum a at time t is 



R ja (t) = H-'a-^t). 



(35) 



It will be convenient to use the scale factor time, /3 = 0. 
Then a(t) = e* and 



dN ja (t) _ 3 3 u 



dt 



An 



(36) 



The comoving volume fractions fj (t) can be found from 
Eq. H31|l . which can be written in a matrix form, 



where f(t) = {fj(t)} and 



(37) 



(38) 



The asymptotic solution of (|37|) at large t has the form 
f(f)=f(°)+se- g * + ... (39) 

Here, f(°) is a constant vector which has nonzero compo- 
nents only in terminal vacua, 



fL 0) - o. 



(40) 



It is clear from Eq. (|33|l that any such vector is an eigen- 
vector of the matrix M with zero eigenvalue, 



Mf = 0. 



(41) 



As shown in Appendix B, all other eigenvalues of M have 
a negative real part, so the solution approaches a con- 
stant at late times. We have denoted by —q the eigen- 
value with the smallest (by magnitude) negative real part 
and by s the corresponding eigenvector. It will also be 
shown in Appendix B that this eigenvalue is real and 
nondenegerate. This follows under the assumption that 
the set of inflating (nonterminal) vacua is "irreducible," 
i.e. it cannot be split into groups of noninteracting vacua, 
where each vacuum from one group never nucleates any 
vacua from other groups. If there were such noninteract- 
ing groups, the present considerations would still apply 
to each group separately. 

The asymptotic values of the terminal components fm 
depend on the choice of initial conditions. For any phys- 
ical choice, we should have fm > 0. Moreover, since f m 



can only increase with time, we must have 
s m < 0. 

At the same time, we should have 



>0, 



(42) 



(43) 



because f a (t) > and f^ = 0. (Here, Eqs. |@2J and 
(|43|l are justified on physical grounds. In Appendix B, 
these equations will be derived more rigorously.) 

It follows from Eqs. (O, and gUJ) that f a w 
Substituting this into H36fl and integrating over t, we ob- 
tain 



-qt 



N ja (t) 



3 rrl 1 

Air a 3-q 



(3-q)t 



(44) 
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The next step is to impose the cutoff. Our prescription 
is to include only bubbles of comoving size greater than 
e. For bubbles of type j nucleating in vacuum a, this 
means that only bubbles nucleated prior to 



t 



(0 

ja 



ln(eH a ) 



(45) 



should be included. With this cutoff, Eq. (|44|l gives 



(eh _ 3 1 



c-(3-3) 



(46) 

As expected, Nj — » oo as e — > 0. Our proposal is that 
Pj oc Nj] hence, 



(47) 



The problem of calculating ^ has thus been reduced to 
finding the dominant eigenvalue q and the corresponding 
eigenvector s. 

We emphasize that our cutoff procedure is indepen- 
dent of coordinate transformations at future infinity. Any 
smooth change of coordinates will locally be seen as a lin- 
ear transformation, which amounts to a constant rescal- 
ing of bubble sizes, accompanied by a rotation. Rescal- 
ing is generally different in different directions, so the 
shapes of the bubbles will be distorted from (approxi- 
mate) spheres to ellipsoids. However, in a sufficiently 
small neighborhood of comoving size S, all types of bub- 
bles are distorted in the same way, so the bubble counting 
should not be affected by the distortion. (For a non- 
spherical bubble, the bubble size can be defined as its 
maximum extent, or the major axis of the ellipsoid.) The 
bubble density is dominated by the bubbles which formed 
at very late times, and which therefore have a very small 
comoving size. Hence, in any neighborhood of size S the 
relative numbers of bubbles will not be affected by con- 
stant rescalings in the limit e — > 0. This means that 
calculation of pj in any such neighborhood will give the 
same result, and it is clear that this result will also hold 
in a comoving region of any size. 



B. Models with quantum diffusion 

In quantum diffusion models, the calculation can be 
performed in a similar way, except the role of bubbles is 
now played by newly thermalized /i-regions. The rate at 
which the comoving volume thcrmalizes can be calculated 
by using the FP equation. 

The FP equation for the distribution function F(<p, t) 
applies in the region of </)-space outside the thermalization 
boundaries si^. The boundary condition on F requires 
that diffusion vanish on si^ In the notation of Sec. 
II.B, 



n a ^-[D{<f)F] 

0<pa 



((f>esi j) ) 



(48) 



(7) 

where n a is the normal to S* . 

The fraction of comoving volume which ends up in a 
thermalized region of type j per unit time with the fields 
(p a * in an infinitesimal surface clement dSi of the ther- 
malization boundary si is given by 



h a (4>*)Ja{<P*,t)dS 



U) 



(49) 



Using Eqs. (|12|l and l|14|) and the boundary condition 
(|48|l . this can be expressed as 



1 

4tt 



H^i^Ficp^dSi^, (50) 



where (dH/dn) = h a (dH/d(p a ). 

The fraction of comoving volume that thermalizes dur- 
ing the scale factor time interval dt ((3 = 0) is thus given 
by 



1 

47r' 



-R- 







dH 
dn 



F{<p*,t)dS*dt. 



(51) 



The asymptotic comoving size of a pocket universe rj is 
much bigger than the comoving size of an /i-region at the 
time of thermalization. Let us denote by A* the ratio of 
these two sizes, 



(52) 



where i* is the time of thermalization. In fact, rj is of the 
order of the comoving size of the horizon at the time t q 
when slow roll begins, which is of order H q e~ tq . Thus, 5 



A* 



(H./H q )e<t' 



-*«) 



(H*/H q )Z, 



(53) 



where Z is approximately given by Q. The number of 
pocket universes which begin thermalizing at time t then 
satisfies 



dNj(<p*,t) 



3ff 2 



.,3/ 



(4tt) 2 A2 



dH 
on 



F(<p*,t)dS*dt. 



(54) 

The asymptotic solution of the FP equation has the 
form 



F(<M) «Fi(0)exp(-git), 



(55) 



where — q\ is the smallest (by magnitude) eigenvalue of 
the FP operator. Substituting this in (|54(l and integrat- 
ing over t, we find the number of /i-regions that have 
thermalized up to time t with <p* in a given surface el- 
ement dS*. The probabilities pj are then found after a 



Regions in the vicinity of an h region which has just entered slow 
roll, have field values which are also close to the boundary where 
slow roll begins, and so they are likely to cross it soon after. Be- 
cause of such correlations, we expect that rj is somewhat larger 
than the comoving horizon size at the time t q , by a factor of a 
few, which can be determined from numerical simulations. 
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cutoff at t {t 



ln[eff ((/>*) /A*] and integration over 



Pj 



Ftfa). (56) 



This equation allows us to calculate pj once we find the 
dominant eigenvalue q\ and the corresponding eigenfunc- 
tion F\{<j>). It should be noted that in the slow roll re- 
gion, the diffusion term can be neglected and the Fokker- 
Planck equation becomes first order. One then readily 
finds that oc Z q (<{>), where Z is the slow roll expan- 
sion factor. With A* given by l|53[) . it is clear that the 
Pj do not depend on where exactly we choose to define 
the thermalization boundary. The integrand can be eval- 
uated anywhere in the slow roll regime, as long as dif- 
fusion is negligible. In fact, Eq. (|56|l is approximately 
valid if we substitute the subindex *, corresponding to 
thermalization, by the subindex q, indicating the onset 
of the slow roll regime. In this case X q ~ 1. Note that 
the value of q\ is typically very small. 

In the preceding subsection we assumed that the 
physics of eternal inflation is described by bubble nu- 
cleation, and in the present subsection we assumed that 
it is described by quantum diffusion. The general case, 
when both mechanisms are present, can be described by 
combining the formalisms outlined in these subsections 

m. 



IV. COMOVING PROBABILITY 

An alternative weight factor for the different types of 
pockets, p^, was used in Ref. 0. This is the so-called 
comoving probability, which can be defined as the prob- 
ability for a comoving observer, starting near the top of 
the potential V(<f>, X), to end up in a pocket of type j. 
Let us now briefly describe, for comparison, the calcula- 
tion of such probabilities. 



A. Models with quantum diffusion 

The form of the FP equation depends on the time 
parametrization parameter /?, and the resulting distri- 
bution F(cf>,t) is, of course, also /3-dependent. However, 
as shown in [13, this distribution can be used to define 
reparametrization- invariant probabilities . 

The probability for a comoving observer to end up in 
a thermalized region of type j is given by 



Pi 



dS*n a ( 



dtJ a (<t>*,t), 



where the first integration is over </>„ £ gam, using 

Eqs. l(T2|) and 1(11)1 and the boundary condition <(lg|) . this 
can be expressed as 



(57) 



1 

4tt 



Pi = -T dS 



?/#*)> 



(58) 



where we have introduced 

V#) = HP- 1 ^) I dtF(<f>,t). 



(59) 



Integrating the FP equation ^li p .(|12 p and the boundary 
condition (|48|l over time, we obtain the corresponding 
equations for ip((f>): 







1 d 



8tt 2 



1 dH 



An < 



-1> 



-F o (0), (60) 



On 



(H 3 iJj) = 



(61) 



Here, Fo(<f>) — F(<f>, 0) is the initial distribution at t = 0. 

The function ip{<j)) is uniquely determined by Eq. I|6U|) 
with the boundary conditions 1)61(1. It can then be used 
in Eq. I|58|) to evaluate the probabilities p^. Note that 



the parameter (3 has been absorbed in the definition (|5t 
of ip((f>) and does not appear in Eqs. (JHEJl , (EOJ > EJ ■ This 
shows that this definition of Pj is independent of time 
parametrization. 

The physical origin of this reparametrization invari- 
ance is easy to understand. If we start with a large in- 
flating volume at t = 0, different parts of this volume will 
thermalize into different minima of the potential, and at 
any time t there will be parts of this volume that are still 
inflating. However, in the limit t — > 00 all the comov- 
ing volume will be thermalized, except a part of measure 
zero. Different choices of time parametrization affect the 
division of the volume into thermalized and still inflating 
regions, but this has no effect on p^, since the inflating 
regions represent only an infinitesimal part of the co- 
moving volume at t — > 00. In other words, the fraction 
of comoving volume that ends up in a given minimum of 
the potential is gauge-independent. 



B. Models with bubble nucleation 

Likewise, in models with bubble nucleation it is conve- 
nient to introduce the new variable 

^ = Y H J~ 4 J Q ( 62 ) 
Integrating Eq. (|31|l over time, we obtain 

/j(oo) - fj(0) = ^(-Xkjipj + Xjktpk) (63) 



In the asymptotic future, all the comoving volume will 
be in the terminal vacua. Hence, 



/a (OO) =0 



(64) 



for all non-terminal vacua, whereas the probabilities for 
terminal vacua are given by 



Pn = /n(«)). 



(65) 
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On the other hand, the initial distribution is concentrated 
at high-energy vacua; hence 

f n (0) = 0. (66) 
Now, it follows from the above equations that 

y^(Aq 7 y> 7 - A 7Q Vq) = / 7 (0) (67) 

a 

and 

Pn =^2K a 1p a - (68) 

a 

The quantities ipa can be determined from (|67|l . and the 
probabilities p° can then be found from Eq. (|68(l . As 
before, the time reparametrization parameter /3 has been 
absorbed in the definition of ip a , and the probabilities do 
not depend on the choice of the time variable. 

C. Problems with Pj 

The probabilities p^ in and lj6l3|l are independent of 
time parametrization, but they do depend on the initial 
probability distribution Fa(<j)) or / Q (0). We assumed that 
the initial distribution is concentrated near the maximum 
of the potential, or in the highest-energy false vacuum. 
However, if the potential has several peaks of comparable 
height, different values of p c - will be obtained starting 
from different peaks. 

It is possible that the initial distribution is to be found 
from the wave function of the universe, which determines 
the probability distribution for the initial states of the 
universe as it nucleates out of nothing. It is well known 
that, although inflating spacetimes are generically eternal 
to the future, all past-directed geodesies in such space- 
times are incomplete, except perhaps a set of measure 
zero (see [24| and references therein) . This indicates that 
the inflating region of spacetime has a boundary in the 
past, and some new physics (other than inflation) is nec- 
essary to determine the initial conditions at that bound- 
ary. The prime candidate for the theory of cosmic initial 
conditions is quantum cosmology, which suggests that 
the universe starts as a small, closed 3-geometry and im- 
mediately enters the regime of eternal inflation. 

If, for example, we adopt the tunneling wave function 
of the universe, the initial distribution is given by |25ll26| 

W)ocexp(-^). (69) 

This distribution favors large values of V((f>). So, if the 
potential is dominated by a single peak, the initial dis- 
tribution will be concentrated at that peak. 

The distribution 1)69(1 may be a plausible choice for 
the initial state of the universe, but the conclusion that 
the probabilities of different observations in an eternally 



inflating universe have some dependence on this choice 
appears to be counter-intuitive. Even though inflation 
must have had a beginning, one expects that once it be- 
gan, the universe will quickly forget its initial conditions. 
Another problem is that, from (|64|l . 

Pa = 0, (70) 

for all non-terminal vacua. If we use for determining 
the probabilities of being in different pockets, all non- 
terminal vacua are given zero probability. However, there 
seems to be no good reason to discard all such vacua. 
Note that the relative abundance of the corresponding 
pockets, discussed in the previous section, can be sizable. 
Also, our own low energy vacuum may have a small pos- 
itive cosmological constant, in which case it would be 
non-terminal. This seems to disfavor the use of p? as a 
relevant weight factor. 

The weight factors Pj which we discussed in Section 
III do not suffer from these problems, and therefore seem 
more suitable for our present purposes. 

V. SOME EXAMPLES 

We shall now use some examples to illustrate the simi- 
larities and differences between the weight factors pj and 
p C y Consider first a very simple model with one false vac- 
uum F and two terminal vacua A and B [FAB model). 
The allowed transitions in this model are shown by the 
"schematic" 

A^F^B. (71) 

The vector f then has three components, f = 
(/f,/a,/b), and the evolution equations i|31|) have the 
form 

= ~{kaf +Kbf)If, (72) 

at 

^ = n AF f F , (73) 

~dt =KBFfF - (?4) 
The first of these equations gives 

f F (t) = Ce~«\ (75) 

with C — const and 

q = k af + kbf- (76) 

With the initial conditions f (0) = {1, 0, 0}, we have C = 
1, and the other two equations yield 

f j (t) = ( KjF /q){l-e-« t ), (77) 

where j = A, B. The comoving probabilities are defined 
by p C j cx fj(oo), which gives 

p) oc Kj F . (78) 
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This is in accord with the intuitive expectation: the prob- 
ability of a vacuum is proportional to the nucleation rate 
of the corresponding bubbles. 

We now compare this with the CHC (comoving hori- 
zon cutoff) approach. The transition matrix M corre- 
sponding to Eqs. I|72 |l -I|74 |l has two zero eigenvalues with 
eigenvectors (0,1,0) and (0,0,1), and a negative eigen- 
value —q [where q is given by Eq. (|76[) ] and eigenvector 
s = (1, —KAF/q, —kbf/q)- For °ur FAB model, the sum 
in Eq. I)47|) has only one term (the one with a = F), and 
the only j-dependent factor in (|47Jl is KjF- Thus, the 
CHC prescription gives 



ijF, 



(79) 



in agreement with Eq. (|78(1 . 

We next consider a model with an intermediate vac- 
uum / between F and B, 



A^F 



B 



(80) 



(FAB I model). Bubbles of A and I nucleate in F, but 
bubbles of I are themselves sites of eternal inflation, with 
an infinite number of B-bubbles nucleating in each of 
them. This model illustrates the difference between the 
weight factors Pj and p^. 

The evolution equations for the FABI model are 



dfF 
dt 
dh 
dt 



-(kaf + K/f)/F) 
-KBlfl + KlFfF, 

kaf If, 



dU _ 
dt 
dfs 
dt 



(81) 
(82) 
(83) 
(84) 



Once again, we start with the comoving probabilities. 
The solution of Eqs. I|81|) - (|84|l with the initial conditions 
(/F,/j,/A,/n)(t = 0) = (l,0,0,0) is 



fp(t) 



-'ii 



f I {t)=C(e-i t -e- KBlt ) 
f A (t) = ^ (1 - e"*) 



f B (t) = 



KlF 



C e 



kbi 



-qi 



where 



and 



q = k A f + k if 



C 



K IF 



kbi - q 

The comoving probabilities for this model are 
/a(oo) kaf 



Pa 
Pb 



Jb(oo) 



k if 



(85) 
(86) 

(87) 
(88) 



(89) 
(90) 

(91) 



This result is easy to understand. We start with all 
comoving volume in F. This volume is then divided be- 
tween the bubbles of A and / in the ratio kaf/^if- All 
the comoving volume in /-bubbles is eventually turned 
into i?-bubbles, and thus the ratio of comoving volumes 
in A and B is given by (|91|) . 

Let us now compare this with the CHC method. The 
matrix M for the FABI model has two nonzero eigen- 
values: — q with q from l|89|) and — q' with q' — Kbi- The 
corresponding eigenvectors are 



s = (sf,si, SA, sb) 



and 



q{q 



K IF q' 



kaf (q ■ 



K IF q' 



1 



(0,-1,0,1). 



(92) 



(93) 



The asymptotic behavior of the model depends on the 
relative magnitude of q and q' . 
For q < q' , or 



kaf + kif < kbi, 
q is the dominant eigenvalue, and Eq. I|47|l gives 



PA 
PB 



H F \ q K AF (q' - q) 



KiFq' 



(94) 



(95) 



Conversely, if q > q', then q' is the dominant eigenvalue, 
and we find 



Pb 



(96) 



Once again, these results are easy to understand. For 
q > q' , the high-energy vacuum F decays faster than the 
intermediate vacuum /. So, in the asymptotic regime the 
comoving volume is divided between bubbles of A and /, 
and each bubble of / still continues to produce bubbles of 
B. As a result, A-bubbles are completely outnumbered 
by £>-bubblcs. In the opposite limit, when q <C q' ~ 1, 
/ turns into B very quickly, and in this sense the model 
is not very different from the FAB model. In this limit, 
Eq. (|95|l gives pa/pb ~ kaf/kif, in agreement with 

ED- 

It is not difficult to extend this analysis to the case 
where the intermediate vacuum / is recyclable, so the 
schematic is A ^ F ^ I ^ B. This is done in Appendix 
A. 

We thus see that both weight factors pi and p\ give 
results which are in agreement with intuitive expecta- 
tions. A priori, it is not clear which one of these two 
objects is more useful for the purposes of defining prob- 
abilities for the constants of nature. Progress can be 
made by working out their values in a variety of models, 
to see which, if any, of the two can be included in a rea- 
sonable definition of probabilities. As discussed in the 
previous section, it appears that on general grounds the 
Pi are preferable. This is because of their independence 
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of the initial conditions, and because they assign non- 
vanishing probabilities to non-terminal vacua (unlike the 
co-moving probabilities). Also, the example considered 
in this section indicates that in models where there are 
intermediate vacua which are also eternally inflating, the 
probabilities pi seem to better represent the actual dis- 
tribution of terminal pockets A and B. Hence, in what 
follows, we shall concentrate on pi. 



VI. THE FULL DISTRIBUTION 

Let us now address the question of calculating prob- 
abilities for observations in the case when there are dif- 
ferent types of pockets. The full probability distribution 
may be written as 

P obs {X) oc^PjiXitj^iX-.tj). (97) 

3 

where n^ s (X;tj) is defined as the number of observers 
that will evolve per unit comoving volume with specified 
values of the fields X, and tj is a local time variable in 

a pocket of type j. The calculation of i%? (X; tj) will 
not be discussed in this paper, and we shall simply as- 
sume that, given the reference class of observers, this 
can be calculated from first principles. The distribution 
Pj(X;tj) should represent the physical volume fraction 
in pockets of type j, with specified values of the fields 
X. Within a given pocket, the quantities Pj(X;tj) and 
n b s (X] tj) should be calculated at the same time tj, but 
it is not important which moment of time and which time 
variable we choose. The volume grows and the density 
decreases with the expansion, but of course the prod- 
uct does not change. (Note that Ti^ s (X;tj) includes 
all observers — present, past, and future — that will ever 
evolve in a comoving volume having unit size at tj , so its 
t-dependence is simply aJ 3 (X;tj), where a,j(X;t) is the 
scale factor.) 

In Ref. |l(| an attempt was made to compare the vol- 
umes in different pockets. The idea was to quantify 
the amount of expansion measured from the inflaton's 
last visit to the highest point in its potential. For that 
purpose, some fictitious markers were introduced, which 
were produced at a constant rate in regions where the 
field is at the top of the potential, and which where sub- 
sequently diluted. The mean separation of the markers 
in the different thermalized regions would then indicate 
how much expansion intervened since the field was at the 
top. This method runs into some conceptual problems. 
If there are several peaks of comparable height in the 
potential, then it is not clear why we should only count 
one of them as the source of markers. Also, since mark- 
ers are only produced at a very special field value, all 
field values which are distant from it get rewarded with 
a large expansion factor, for no particularly compelling 
reason. In the face of such difficulties, one option would 



be to content ourselves with the volume distributions for 
individual pockets as the last frontier of predictivity. 

Nevertheless, the structure of the eternally inflating 
spacetime is well understood, and we feel that a useful 
characterization of the relative likelihood the constants of 
nature in pockets of different types can be obtained from 
the distributions pj and P(X;j) which we have consid- 
ered in the previous sections. 

Let us start, for simplicity, with the case where the 
different pockets are generated by bubble nucleation. 
Suppose also that different pockets have different val- 
ues of the constants, but that these constants do not 
vary within a given bubble. Our goal is to determine the 
probability of being in a given pocket j, P^/ s - In this 
case, the problem is that of combining pj and into 

Poll [f° r t ne time being, there is no P(X;j) to worry 
about, since the constants X do not vary within a given 
pocket] . 

For that, we need to define a co-moving reference scale 
specifying the size of the regions in which the observers 
are to be counted, and we can write 

Pl^PiR^lit), (98) 

where Rj(t) is the physical size of the reference scale 
and we have omitted the subscript j of tj. To deter- 
mine Rj{t), we note that at early times, the dynamics of 
the open FRW universe inside a bubble is dominated by 
curvature, and we have a,j ~ t for all types of bubbles. 
Later on, when the curvature scale grows to the size of 
the Hubble radius associated to the corresponding local 
energy density, the curvature dominated phase gives way 
to slow roll inflation inside the bubble. For a quasi-de 
Sitter slow roll phase we would have 

aj{t)^{l/H^) S mh{H^t), (99) 

where H q is the Hubble rate at the beginning of the slow 
roll. The specific form of the scale factor at late times 
is not important for our argument. The point is that 
for times much smaller than the onset of slow roll, t <C 
1/Hq , all pockets are nearly identical, with scale factor 
dj{t) « t. This suggests that the reference scales should 
be chosen so that Rj (t) is the same in all pockets at some 

t = e < l/H^\ Then, up to a constant, 

Rj(t)= aj (t). (100) 

For times t ^> 1/Hq this can be expressed as Rj(t) « 

(l/Hq^)Zj(t), where Zj is the redshift factor from the 
time at which slow roll begins. Alternatively, we note 
that Rj(t) is the comoving curvature scale, which can be 
defined without any reference to E g . 

Now, this can be straightforwardly generalized to the 
case when there are some continuous variables X, by sim- 
ply including the factor P e (X;j) — the normalized distri- 
bution for X on the surface t = e: 

P obs (X) cx Y,P3HXj)R Z J {X)7i^ s (X). (101) 

3 
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The volume distribution of X on the surface E g is 



P q (X;j) 



L P £ (X;j)H^(X), 



(102) 



where we abreviate Hj = H q , and is the normaliza- 
tion factor 

^ = / P t {X;j)Hf{X)dS q = (Hf) e . (103) 

Combining the above equations, the fraction of volume 
in pockets of type j with values of the constants X takes 
the form 



/',i.Vi xp jVj P q (X;j)Z?(X). 



(104) 



Here, P g is the single pocket volume distribution nor- 
malized at the onset of slow roll (in the case of diffusion, 
this can be estimated as the ergodic factor P q (X;j) ~ 
N exp[S(X; j)], or using Eq. 124J1 in the case of bubble 
nucleation). The factor Vj has the dimension of volume. 
From (jTU2|) we have VjH^P q = P e . The left hand side 
of this equation is normalized to unity, and therefore we 
may write Vj in terms of the volume distribution on S g , 



Hf(X)P q (X;j)dS q 



(105) 



Note that the product P h (X;j) oc H^(X)P q (X;j) featur- 
ing in the denominator of Vj is the number distribution 
of /j-regions on E g with the given value of X in pockets 
j. An alternative form of (|104fl is given in terms of this 
number distribution by 



Pj(X) <x Pj P h {X;j)Hf{X)Z?{X), 



(106) 



where Ph is normalized on S q . This result has a simple 
intuitive interpretation. The volume distribution is pro- 
portional to the relative number of pockets of a given type 
Pj, times the relative number of Hubble regions within 
these pockets with given values of X, Ph(X;j), on the 
surface where slow roll starts, times the volume of these 
regions H~ 3 , times the subsequent growth factor Z|. 

Our derivation of (|104fl . or its equivalent form (|106fl . 
has been based on the case of bubbles, but it is possi- 
ble to adopt the same prescription also for the case of 
pocket universes generated by diffusion, since the result 
only makes reference to the surface £ q where slow roll 
begins and to the subsequent expansion thereof. 

At present, we cannot claim that the prescription we 
have introduced here is the only consistent one. How- 
ever, the hope is that from the careful analysis of a suf- 
ficiently wide range of possibilities a unique prescription 
will emerge. To illustrate this point, let us consider two 
alternatives which on the surface may seem reasonable. 

Let us first consider a possibility which we shall refer 
to as Alternative A. This consists of identifying pj with 
the volume fraction in pockets of type j. In this case, we 
would have 



Here, P(X;j) is the distribution P(X;j) normalized on 
the corresponding thcrmalization boundary, 



dS^P(X;j) 



(108) 



Although the choice (|107H may appear reasonable at first 
sight, it is unsuitable for substitution into Eq. I|97|) . In- 
deed, the product Pj(X)n ! obs {X) should be independent 
of the time tj at which we choose to calculate it. But with 
the present choice, Pj(X) stays constant, while n 3 obs (X) 
dilutes with the expansion. 

To avoid this problem, we may then think of another 
possibility, which we shall refer to as Alternative B, 
where 



p obs (x) = J2pj p L(x)- 



(109) 



Here, it is Pl bs {X) = NP{X]j)n' obs which is normal- 
ized to one. With this, the time at which we calculate 
P(X: j)n J obs is unimportant, since this product is inde- 
pendent of time. This proposal, however, does not take 
the form (|97() and contradicts our starting point that all 
members of the reference class of observers carry equal 
weight (which is in fact the very definition of reference 
class). Suppose we have a single false vacuum decaying 
into two types of bubbles with equal nucleation rates, 
so that pi = P2- The vacuum in type-1 bubbles is very 
hospitable to observers, while the vacuum in type-2 bub- 
bles is almost lethal and density of observers is strongly 
suppressed. One expects that observers are much more 
likely to find themselves in a type-1 bubble. However, af- 
ter the normalization prescribed in l|109|) , we would con- 
clude that we are equally likely to be in either type of 
bubble. The basic problem here is that the factors pj are 
determined by the transition rates, or by the diffusion 
dynamics, and therefore they will be insensitive to the 
amount of slow roll inside of a given pocket — and much 
less to the number of observers which will subsequently 
develop. Hence, we must also reject this alternative. 

The preceeding discussion illustrates that, as we have 
seen in the case of Alternatives A and B, certain 
reasonably-looking prescriptions can be discarded after 
careful consideration of their consistency. Another way 
of narrowing down the possibilities is by working out 
their predictions and comparing them with the data. 6 Al- 
though a full analysis is beyond the scope of the present 



P j (X)=p j P(X;j). 



(107) 



6 This type of analysis has already proved useful when applied to 
definitions of probability based on a constant-time cutoff. For a 
wide class of time variables, such definitions lead to absurd pre- 
dictions that we should find our selves in a deep, nearly spherical 
gravitational potential well ll'Ttl and that the CMB tempera tur e 
should be many orders of magnitude higher than observed Il4l . 
(As we mentioned in the Introduction, another objection against 
constant-time cutoffs is that they rely on an arbitrary choice of 
the time variable.) 
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paper, we believe that the definition i|104|) [or equiva- 
lently 1)11) is the most appealing amongst those which 
we have considered, and has the best chance of being 
correct. 



VII. SUMMARY AND DISCUSSION 

In the context of theories with many vacua, such as 
the landscape of string theory, the low energy constants 
of Nature are randomized during inflation, and will vary 
from place to place in the thcrmalized post-inflationary 
regions. An interesting question is then to find the vol- 
ume distribution of the constants in such thermalized 
regions. In this paper, we have developed some tools 
towards this goal. 

The problem splits into two parts. First of all, we 
should find the distribution of constants within a given 
pocket universe. Following up on the approach of Refs. 
H, H|, we have motivated an analytic estimate for this 
distribution. According to Eq. J2J, the internal volume 
distribution in a given pocket of type j takes the form 

P(I;j)c(P,(I;j)Z 3 (I;j), (110) 

where P q (X;j) and Z(X;j) are the distribution at the 
onset of slow roll and the slow roll expansion factor re- 
spectively. We argued in Section II. C that 

P q (X;j) « H- 2 (X;j)exp[S(X;j)]. (Ill) 

This is what we called the ergodic conjecture. 7 It should 
be understood that our justification of (|lllf) has been 
rather heuristic, and the limits of its validity should be 
further explored. 

Next, we have introduced a weight factor pj which 
counts the relative number of pockets of type j. The 
counting is done at the future boundary of spacetime, 
by considering all pockets of comoving size larger than 
some e, and then letting e — > 0. We have called this the 
comoving horizon cutoff (CHC), since the comoving size 
of a pocket is roughly given by the size of the horizon at 
the time when the pocket is formed. As we explained at 
the end of Section III. A, such a cutoff is independent of 
coordinate transformations at the future boundary. 

For comparison, we have also discussed an alternative 
weight factor p^, defined in 01 as the fraction of comov- 
ing volume which ends up in pockets of type j. Unlike 
the pj's discussed above, the comoving volume fraction 
p^ strongly depends on initial conditions. Moreover, all 
co-moving volume ends up in terminal vacua (these are 
vacua which cannot further decay into other vacua) . Be- 
cause of that, all non-terminal vacua are assigned zero 



7 In models of bubble nucleation followed by slow roll inside the 
bubble, Eq. 11111 should be replaced with an expression of the 
form I2H . 



probability. However, there seems to be no good reason 
to discard all such vacua, since the relative abundance of 
the corresponding pockets can be sizable. Also, our own 
low energy vacuum may have a small positive cosmolog- 
ical constant, in which case it would be non-terminal. 
Because of these two problems, the probabilities p C j seem 
less relevant for our present purposes. 

In the case where all terminal vacua emanate directly 
from a single eternally inflating false vacuum, we find 
that the weight factors pj and p^ agree with each other, 
and they are basically proportional to the correspond- 
ing nucleation rates. The differences between pj and p^ 
are illustrated in a model where there is an intermediate 
false vacuum which is also eternally inflating. Again, the 
weight factors pj seem to give a better representation of 
the actual distribution of pockets in the multiverse. 

We have combined the above distributions into an ob- 
ject describing the thcrmalized volume distribution of the 
constants of nature in the inflationary multiverse with 
pockets of different types. If these constants are deter- 
mined by the values of certain light fields X, then we 
have argued that the volume fraction occupied by values 
X in pockets of type j is given by 

dPj(X) = Pj(X) dS& = Pj dS®. 

(112) 

Here, Sq is the boundary (in field space) where quan- 
tum diffusion turns into classical slow roll. The factor pj , 
whose calculation is discussed in Section III, takes into 
account the relative numbers of pockets of type j. The 
quantity Ph(X;j) is the number distribution of Hubble- 
size regions with values of the fields X within pockets of 
type j at S 9 (see Fig. 2 for the space-time structure of a 
pocket universe), and H~ 3 (X) is the volume of one such 
region. According to the ergodic conjecture, this distri- 
bution is essentially given by Eq. (|lllfl [or Eq. (|24ll in 
the case of bubbles] . The hat on Ph indicates that it is 

normalized as / P h (X;j)dS q 3) = 1. Zf(X) is the expan- 
sion factor from the point X on S q to the thermalization 
boundary S*, along the classical slow roll trajectory. 

Although our proposal for comparing the probabilities 
in pockets of different types seems to be well motivated, it 
may not be the unique possibility. Settling this issue may 
require further research. For comparison, we have also 
explored some alternatives. In the first one, the weight 
factor pj is formally identified with the volume fraction 
in pockets of type j, and in the second one, pj is formally 
identified with the fraction of observers in pockets of type 
j. We have shown that both of these options would lead 
to inconsistencies. 

With the analytic approximations we have suggested, 
the distribution (|112f> can be readily calculated in a vari- 
ety of models. The probability distribution for the mea- 
surements of observers in a given reference class can be 
obtained by substituting i|112[l into Eq. I|97|) . Application 
to examples which may be relevant to the landscape of 
string theory (such as the Bousso-Polchinski scenario 1] ) 
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is left for further research. 
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Note added 

After this paper was submitted, Easther, Lim, and 
Martin (ELM) suggested an alternative prescription for 
computing the weights pj assigned to different types of 
pocket universes |3(j. In this note we show that their 
prescription is equivalent to the CHC method. 

The ELM prescription is formulated for a model with 
bubble nucleation, and we shall restrict our attention to 
that case. The main problem of counting bubbles is to 
select a large but finite subset from the infinite set of 
all the bubbles created during the infinitely long evolu- 
tion. Once a finite subset is selected, the probability ratio 
Pj/pk is found as the ratio of the number of bubbles of 
types j and k from the subset, in the limit of large sub- 
set size. The ELM proposal consists of choosing a large 
but finite number of randomly drawn comoving world- 
lines and selecting the subset of bubbles that intersect at 
least one of these worldlines. The distribution of world- 
lines is assumed to be statistically independent of the 
bubble nucleation process. Then the resulting weights pj 
are independent of the choice of worldlines. 

The subset of bubbles selected by the ELM prescrip- 
tion differs from the subset selected by the CHC method 
because some bubbles in the ELM subset may have an 
arbitrarily small comoving size. However, we shall now 
show that the ELM prescription produces the same re- 
sult for Pj/pk as the CHC method, for a generic bubble 
nucleation model (with or without recycling). 

The above mentioned assumption of statistical inde- 
pendence ("the worldlines do not know about the pock- 
ets" ) is equivalent to assuming that there exists a well- 
defined probability density R(j;t)dt for a randomly cho- 
sen worldline from the congruence to end in a bubble of 
type j at a time t. Let us draw the worldlines uniformly 
with a constant probability per unit comoving 3- volume; 
this is an admissible distribution in the sense of the ELM 
prescription. Then R{j\t) is equal to the probability of 
nucleating a bubble of type j at time t at a randomly 
chosen comoving point. It is clear from Eq. i|31|) that 

R(j;t) = J2*iMt), (113) 

a 

where the index a spans only the non-terminal vacua. 

To simplify the calculation, we consider a uniform cu- 
bic grid of worldlines where all the cubes have sides e / \/3. 
Then the smallest and the largest (comoving) distances 



between neighboring worldlines are and e. This is 

again an admissible congruence of worldlines, chosen in- 
dependently of the bubble nucleation process. The subset 
of bubbles selected by the ELM prescription will include 
all bubbles of comoving size L > e (the counting of which 
constitutes the CHC prescription), as well as some bub- 
bles of smaller comoving size. So it remains to show that 
the additional counting of bubbles of smaller size will not 
influence the CHC result, in the limit e — > 0. 

The total number of selected bubbles with sizes L < e 
is no greater than the number of worldlines that do not 
intersect larger bubbles of size L > e. It follows from 
Eqs. (|113[) . and the condition /" = 0, that the prob- 
ability R(j; t) of a worldline encountering a bubble at late 
times t is exponentially small, cx e~ qt . Due to Eq. 145|l . 
the limit e — > corresponds to the limit t — ► oo. Hence, 
the additional counting of smaller-size bubbles produces 
an exponentially small correction to the CHC counting, 
and this correction disappears as t — > oo. Therefore, we 
have recovered the CHC bubble counting result in the 
limit e — > 0. Since the ELM prescription is insensitive to 
the distribution of the worldlines, we conclude that the 
two prescriptions always yield identical results. 



Appendix A: FABI model with recycling 

We now consider a version of the FABI model which 
includes recycling between vacua F and / as indicated in 
the schematic A *— F <-> I — > B. 

Writing out the master equation, we obtain: 



dfF 
dt 



{kaf - 


\- k/f)/f - 


V KFlfl 


(114) 


(kfi - 


I" K B l)fl "I 


- kifJf 


(115) 




df A _ 
dt 


hafJf 


(116) 




dfo _ 
dt 




(117) 



We will first consider the comoving probabilities. Not- 
ing that we take n mn to be independent of time, we can 
integrate equations l|114|) to l|117|) from some initial time 
t = to t = oo to obtain 

- /f(£ = 0) = -(kaf + kif)^f + kfi4>i, (H8) 
-fi(t = 0) = -(kfi + KBl)i>I + kif"4>f, (119) 

f A (t = oo) = KAFtpF, (120) 

f B (t = oo) = Kbi^i, (121) 

where we have defined 

/>oo 

A = / fidt. (122) 
Jo 
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For simplicity we also define 



Using Eq. iBTjl we find as c — > 0, 



kaf = a 
Kjp = b 

Kpj = c 

K B i = d 



(123) 
(124) 
(125) 
(126) 



Solving this system of linear equations for ipi in terms 
of the initial conditions and the transition rates, we find 
the ratio of comoving probabilities as 



p^ = f A {t = ™) = a(c + df F {t = 0)) 
p% f B (t = oo) d(b + af I {t = 0)) , 



(127) 



with a,b,c,d related to the transition rates as given in 

Eqs. (jnsi-inni- 

For the initial condition (fF,fi,fA,fs)(t — 0) = 
(1,0, 0, 0), we have 



p c A _ a(c + d) 



db 



If we consider c <C d, we find 

Pa_ ^ a _ kaf 
P C B b k if 



(128) 



(129) 



in agreement with Eq. (|91[l. 

We now compare this result with the CHC method. We 
solve for the eigenvalues and eigenvectors of M, yielding 
explicit expressions for the Sj. The result is 



s F 



sa 



_ 1 
~ 2 

-q 

2bd 



(a + b 



— ( 

2bd y 



c + d- 


n) 


(130) 


c — d — 


K) 


(131) 


si -- 


q 
d 


(132) 


c-d- 


n) 


(133) 


sb = 


-l 


(134) 



where 



K = y/(a + b- c- d) 2 + 4bc, 



(135) 



— q is the dominant eigenvalue, and s = (sf, si, sa, sb) is 
the corresponding eigenvector. For completeness we note 
that the other nonzero eigenvalue (call it — q') is given by 



q = ^(a + b + c + d + TZ) 



(136) 



and that this differs from q only by the sign of the square 
root term. 

If we now want to calculate the probabilities in the 
limit that c <g; a,b,d (that is, we assume that recycling 
from / back into F is nearly negligible), we can Taylor 
expand sa about c = 0. We find 



a 



ca 



SA ~2bd {a + b - d - Via + b - d]2) - d(a + b-dy 



PA 
PB 

if a + b > d and 

Pa r 

PB 



Hf^ 

Hi 



d(a + b) 



Hi 



bd 



(-a-b + d) 



(138) 



(139) 



if a + b < d, in agreement with the results of Section V, 
where c was neglected from the very beginning. 



Appendix B: Nondegeneracy of the subleading 
eigenvalue 

In this appendix we study the mathematical properties 
of the matrix M defined by Eq. I|38|) . Formulations and 
proofs of some technical lemmas will be deferred until 
the end of this appendix. We shall show that all nonzero 
eigenvalues of M have negative real parts, and that the 
subleading eigenvalue — q of the matrix M is real and 
nondegenerate, under the following assumptions: 

I. The set of all the inflating (nonterminal) vacua can- 
not be split into disconnected groups, where each vacuum 
from one group never nucleates any vacua from other 
groups. 

II. Transitions between any two inflating vacua are re- 
versible: If K mn 0, then K nm 7^ as well. 

III. There exist transitions to some terminal vacua 
with nonzero rates. 

We begin by considering the vector of comoving volume 
fractions f which satisfies Eq. 1)3 70. 

j t f = Mf , (140) 
where the matrix elements of M are defined by 

Mij = Kij - Sij V, K n- (141) 



Conservation of the comoving volume requires that the 
sum of any column of M be equal to zero. It follows from 
this condition that all eigenvalues of M have nonpositive 
real parts (see Lemma 1 below). 

Suppose that there are rij inflating vacua and nt ter- 
minal vacua. We can split the vector f into a direct sum 
of two vectors, and fv t \, representing the volume frac- 
tions in these two types of vacua. Then Eq. I|140|) can be 
rewritten as a system of two vector equations, 



dt 



^ f (*)- Sf w> 



where the matrix R describes the transition rates be- 
tween the inflating vacua and S describes the transition 
rates from inflating to terminal vacua. Thus the matrix 
M has the following block appearance, 



(137) 



M = 



R 
S 



(142) 
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Note that the matrix R is square with dimensions n,- xn,-, 
while the matrix S has dimensions x n t and may not 
be square. 

We already know that the matrix M has a zero eigen- 
value and no positive eigenvalues. We shall now prove, 
under the assumptions I-III above, that the eigenvalue of 
M with the algebraically largest real part is a nondegen- 
erate and negative eigenvalue — q < 0. 

It follows from the assumptions I and II that any in- 
flating vacuum will eventually nucleate bubbles with any 
other inflating vacuum (either directly or after passing 
through bubbles of other inflating vacua). The tech- 
nical term for this property is that the matrix R is 
irreducible: for any indices i,j there exist a chain of 
indices fci, fc2, k s such that all the matrix elements 
Rik^Rk^, ■■■,Rk s j are nonzero. 

Let us now analyze the eigenvalues of M in general. 
Any eigenvector of M has either all zero inflating com- 
ponents, i.e. it is a vector of the form f = (0,f( t j), or 
f = (f(i),f( t )) with iu\ 7^ 0. All vectors of the form 
(0,f( t )) are obviously eigenvectors of M with eigenvalue 
0, while eigenvectors of the form f = (f(i), f(t)) such that 
Mf = Af must satisfy Rf(i) = Af(j). We shall shortly 
demonstrate that all eigenvalues of the matrix R have 
negative real parts and that the maximal eigenvalue — q 
of R is real, negative, and nondegenerate, while the cor- 
responding eigenvector r can be chosen with all posi- 
tive components. It will then follow that the subleading 
eigenvalue of M is equal to —q and is nondegenerate, and 
the corresponding eigenvector of M is s = (r, — (/ _1 Sr). 
Note that the components of the matrix S are nonneg- 
ative, and thus the inflating components of the vector 
s have the opposite sign as compared with the terminal 
components. This confirms our earlier statements <|42[) 
and gSJl. 

Now we shall show that all eigenvalues of the matrix R 
have strictly negative real parts. Since R is irreducible 
and has nonnegative off-diagonal components, there ex- 
ists a nondegenerate eigenvalue Ao of R with the largest 
real part (the maximal eigenvalue) and all other eigen- 
values have a smaller real part (see Lemma 1). By con- 
struction, the matrix R is of the form (|14H . where the 
indices m,n are restricted to the range 1 < (m, n) < n,-. 
It follows from the assumption III that the sums of each 
column of R are not all equal to zero and thus 



i=i 



rii+nt 



Kii < 0. 



(143) 



i=ni-\-\ 



Then it follows from the estimate l|150|) of Theorem 1 
that Ao < 0. However, we would like to prove the strict 
inequality Ao < 0. To this end, we shall assume that 
Ao = and arrive to a contradiction. If Ao = 0, there 
exists an (n^-dimensional) eigenvector v with all positive 
components (Lemma 1) such that Rv = 0. We can now 
express the matrix R as a difference of a matrix P having 
vanishing sums of each column, and a diagonal matrix Q 



whose entries are nonnegative, namely 

Ru=Pi3-Qij> (144) 

m n r 

I', j = Kij ~ S lj ^ /.';,. ]T Pij = 0, (145) 

1=1 i=l 
Qij = <% E Kli - ( 146 ) 

Here the matrix P describes the nucleation of inflating 
vacua and the matrix Q describes the nucleation of ter- 
minal vacua. The equality Rv = then yields 



Pv = Qv > 0. 



(147) 



Since all the components of v are positive, and since 
Qij > and Q =/= by assumption III, we conclude that 
not all components of the vector Pv are zero: 



Pv f 0. 



(148) 



Note that the matrix P definitely has a zero eigenvalue 
since the sums of all columns vanish. Applying Lemma 1 
to the matrix P, we find that the maximal eigenvalue of P 
is Ap < 0, and since the matrix P has a zero eigenvalue, 
this eigenvalue must be maximal, i.e. Ap = 0. Intuitively, 
we expect that this maximal eigenvalue is diminished 
when we subtract a nonnegative and nonzero matrix Q 
from P. Indeed, Lemma 3 shows that the existence of a 
positive vector v satisfying Pv > means that v is the 
eigenvector of P corresponding to Ap = 0, in other words, 
we have Pv = 0, which contradicts Eq. I|148fl . Therefore 
the matrix R cannot have an eigenvalue Ao = 0. 

Finally, we list some statements used in this appendix. 
A matrix (or a vector) is called nonnegative if all com- 
ponents are nonnegative. We write a > b for vectors if 
a, > bi for all i. 

Theorem 1: A nonnegative matrix j4y has a real 
eigenvalue Ao > such that all other eigenvalues Aj, 1 < 
i < n — 1 are smaller than Ao in magnitude, |A$| < Ao. 
If the matrix Ay is irreducible, then the eigenvalue Ao is 
nondegenerate and the corresponding eigenvector can be 
chosen with all positive components. Furthermore, if we 
denote 

*i E 1 - ( 149 ) 

i 

then the largest eigenvalue Ao is bounded by 

mindj < Ao < maxcrj. (150) 

3 3 

This is the Perron- Frobenius theorem and its corollary. 
For a proof, see (2^, chapter 9, and also [2^, chapter 1. 

Lemma 1: If a matrix Mij is irreducible and My > 
for i =/= j, then there exists a nondegenerate real eigen- 
value Ao of Mij such that the corresponding eigenvector 
has all positive components, and all other eigenvalues 
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have smaller real parts, Re A^ < Ao (maximal eigenvalue). 
Moreover, if £\ M ij < for au h th° n < 0. 

Proof: We can choose a real number r such that the 
auxiliary matrix Ay = My + rJy is nonnegative. The 
matrix Ay is irreducible and thus, by Theorem 1, has 
an eigenvalue Ao such that the corresponding eigenvector 
has all positive components^ and all other eigenvalues Xi 
lie within the circle |Aj| < Ao; hence, ReA^ < Ao. Since 
the eigenvalues of My are A^ = A^ — r, the first statement 
follows. The second statement follows from the estimate 
in Theorem 1, which yields 

Ao < max Ay = r + max My < r, (151) 

i i 

hence Ao = Ao — r < 0. 

Lemma 2: If Ay is a nonnegative irreducible matrix, 



and if Ao is the maximal eigenvalue of Ay from Lemma 1 , 
and if there exists a nonnegative vector v ^ such that 
Av — Aov > 0, then v is an eigenvector of A with eigen- 
value Ao. 

A proof of this technical statement is contained in the 
proof of Theorem 4.1 in Ref. |2{|, chapter 1. 

Lemma 3: If a matrix My is irreducible and My > 
for i j, and if Ao is the maximal eigenvalue of My 
from Lemma 1, and if there exists a nonnegative vector 
v 7^ such that Mv > Aov, then v is an eigenvector 
with eigenvalue Ao- 

Proof: We consider the auxiliary nonnegative matrix 
A = M + rl as in the proof of Lemma 1. It follows that 
Av — (Ao + r) v > 0. We note that Ao + r is the maximal 
eigenvalue of A. Now Lemma 2 can be applied to the 
matrix A and it follows that v ^ is an eigenvector of 
A with eigenvalue Ao + r. 
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